ABSTRACT
Bovine coronavirus (BCoV) is a member of pathogenic Betacoronaviruses that has been circulating for several decades in multiple host species. Given the similarity between BCoV and human coronaviruses, the current study aimed to review the complete genomes of 107 BCoV strains available on the GenBank database, collected between 1983 and 2017 from different countries. The maximum-likelihood based phylogenetic analysis revealed three main BCoV genogroups: GI, GII, and GIII. GI is further divided into nine subgenogroups: GI-a to GI-i. The GI-a to GI-d are restricted to Japan, and GI-e to GI-i to the USA. The evolutionary relationships were also inferred using phylogenetic network analysis, revealing two major distinct networks dominated by viruses identified in the USA and Japan, respectively. The USA strains-dominated Network Cluster includes two sub-branches: France/Germany and Japan/China in addition to the United States, while Japan strains-dominated Network Cluster is limited to Japan. Twelve recombination events were determined, including 11 intragenogroup (GI) and one intergenogroup (GII vs. GI-g). The breakpoints of the recombination events were mainly located in ORF1ab and the spike glycoprotein ORF. Interestingly, 10 of 12 recombination events occurred between Japan strains, one between the USA strains, and one from intercontinental recombination (Japan vs. USA). These findings suggest that geographical characteristics, and population density with closer contact, might significantly impact the BCoV infection and co-infection and boost the emergence of more complex virus lineages.